Microscopic derivation of self-consistent equations of Anderson localization in a 
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We present a microscopic derivation of self-consistent equations of Anderson localization in a 
disordered medium of finite size. The derivation leads to a renormalized, position-dependent diffu- 
sion coefficient. The position dependence of the latter is due to the position dependence of return 
probability in a bounded medium. 



INTRODUCTION 



The phenomenon of Anderson localization [l[ has been studied both experimentally and theoretically for already 
half a century [2, 3,0]. It takes place for waves in strongly disordered media when interference effects become plethoric 
in the multiple scattering process. Theoretical description of Anderson localization reached a decisive stage in the 
eighties with the self-consistent (SC) theory of Vollhardt and Wolfle However, in its original form, this theory 
did not fully account for finite-size effects. Later, Van Tiggelen et al. proposed a natural generalization of SC theory 
^ ' to media of finite size by introducing a position- dependent diffusion coefficient D ,Q\. This generalized SC theory has 
. been recently used to study the dynamics of Anderson localization in quasi-one-dimensional 0] and three-dimensional 
Q systems. Meanwhile, the generalized SC equations of Refs. [H, 0, H have never been derived microscopically. Such a 
derivation is highly desirable for at least two reasons. First, our recent results indicate that the position dependence 
I • of D is crucial for the internal consistency of the theory itself and that some of the important features of Anderson 
localization (like the 1 /L 2 scaling of the transmission coefficient with the size L of a disordered sample at the mobility 
edge) cannot be reproduced without fully taking it into account @ . Second, the very fact that D should be position- 
dependent can be questioned in favor of momentum [loj or time [ll[ dependencies studied in the past, unless the 
position dependence of D is given a microscopic justification. This calls for a rigorous derivation of SC equations in a 
medium of finite size, showing the emergence of position-dependent D from microscopic equations of wave propagation 
and clarifying the physics behind it. 

In this paper we present a derivation of SC equations of localization in a finite medium of size L much exceeding the 
y—{ \ two main "microscopic" length scales of the problem: the wavelength A and the mean free path £ due to disorder. Our 
derivation is based on the "Hikami box" formalism HH • We work in the framework of classical wave scattering, 
but our results can be extended to quantum particles (e.g., an electron or an atom at low temperatures) described 
,_H ■ by Schrodinger equation with a disordered potential. Whereas electronic properties of disordered systems have been 
a subject of intense studies over several decades [l|, [j| E3, El, E2, EH EH El) EE E3) the behavior of coherent atomic 
, ensembles (Bose- Einstein condensates) in disordered optical lattices has come into focus only recently pH El) HI, HH • 
Mathematically, the finitcness of the medium comes into play when we evaluate interference corrections to the sum 
of ladder diagrams. These interference corrections are due to infinite series of maximally-crossed diagrams that we 
insert inside the ladders. In the presence of time-reversal invariance, the final result depends on the probability for 
the wave (or quantum particle) to return back to a given point r. Whereas, due to the translational and rotational 
invariance, this return probability is a position-independent quantity in the infinite medium, it becomes position- 
dependent in a medium of finite size. In an open medium, the return probability decreases when the boundary of the 
medium is approached because of the increased probability for the wave to leave the medium through the boundary. 
This leads to a position-dependent renormalized diffusion coefficient D, the renormalization being less important 
near the boundaries of the disordered medium. The dependence of D on r is not known in advance, but has to be 
determined self-consistently by solving a diffusion equation containing the same D. 

Anderson localization is often defined as an asymptotic property of eigenstates of disordered wave (Schrodinger, 
Helmholtz, etc.) equations. The states are said (exponentially) localized if their intensity decays exponentially at 
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large distances. Another widespread definition of Anderson localization is vanishing of the diffusion coefficient [4|. 
Strictly speaking, none of these definitions can be directly applied in open media of finite size, which were a subject of 
extensive work initiated by Thouless and culminated in the scaling theory of localization fiTfl. We are not intended 
to give a review of this work here and the interested reader can find more details in Refs. @, [Jlj]. For our purposes it 
will be sufficient to think of "Anderson localization in a medium of finite size" as of an interference, wave phenomenon 
that would give rise to "truly localized" states if the medium were extended to infinity. 

The paper is organized as follows. In Sec. [II] we review SC theory of localization in infinite and finite media. The 
main "building block" of our derivation — an "interference loop" that we insert inside ladder diagrams to account for 
interference effects in the intensity Green's function — is calculated in Sec. IIIH In Sec. IIVI we sum an infinite series 
of diagrams for the intensity Green's function and obtain the SC equations of localization. Section [V] is devoted to 
boundary conditions and a discussion of energy conservation. Finally, we summarize our main results and discuss 
their implications in Sec. IVII Technical details of calculations are collected in 4 appendices. 



II. THEORETICAL FRAMEWORK 



We consider propagation of a scalar, monochromatic wave of circular frequency w in a disordered three-dimensional 
medium of finite size. The amplitude Green's function G(r, r',oj) obeys the Hclmholtz equation: 

[A r + fc 2 (l + ,i(r))]G(r, r', w ) = S(r - r'). (1) 

Here /i(r) = Se(r)/e is the relative fluctuation of the dielectric constant e(r) = e + Se(r), e is the average dielectric 
constant, k = ^/icu/c is the wave number, and c is the speed of wave in a homogeneous medium with e — 1 (vacuum). 
We assume that //(r) obeys the white-noise Gaussian statistics: 

fc 4 (Mr)MO> = ya(r-r'), (2) 

where angular brackets denote averaging over realizations of disorder and £ is the scattering mean free path. The 
average amplitude Green's function can be calculated assuming weak disorder (k£ 3> 1) Q: 

<G(r,r») = - i7 ^cx P (ik\r-r'\ - ^f) . (3) 

Although this result has been obtained for the infinite medium, it holds in a medium of finite size as well, provided 
that the points r and r' are at least one mean free path from the boundaries. 
In this paper we will be interested in the average intensity Green's function: 

G(r,r',Q) = — <G(r,r> 1 )G*(r,r', Wa )>, (4) 
c 

where u>i = loq + 0/2, oj 2 = w — 0/2, and we omit the dependence of G on the carrier frequency loq. We assume the 
latter to be fixed in the remainder of the paper. Physically, the Fourier transform C(r, r', t — t') of Eq. (HJ) describes 
the density of wave energy at r at time t due to a short pulse emitted at time t' by a point source at r'. For a 
quantum particle, G can be interpreted as a probability density of finding the particle in the vicinity of point r at 
time t, provided that the particle was at r' at time t' ("probability of quantum diffusion") [l4T ]. 

The analysis of the intensity Green's function is generally complicated and relatively simple results can be obtained 
only for weak disorder (ki ^> 1) at large spatial scales (|r — r'| ^> I) and for slow dynamics (£1 <C cuq, c/£). Under 
these assumptions, one derives the diffusion equation for the intensity Green's function [3I [T3|: 

{-m-D B A r )C(r,r',Q.) = 5(r-r'), (5) 

where Db — c£/3 is the Boltzmann diffusion coefficient. This equation holds in the infinite as well as in finite media, 
provided that it is supplemented with appropriate boundary conditions 0, [13, H3] in the latter case. Obviously, Eq. 
([5]) ignores interference effect and treats the wave as a classical particle that propagates through a disordered medium 
by diffusion. Vollhardt and Wolfle [5[ have shown that interference effects lead to a renormalization of Db in Eq. 
The renormalized diffusion coefficient D(Q) obeys [30j : 

1 , 6tt f dQ 1 

(6) 



D(fl) D B k 2 £ J (2tt) 3 -iO + D(0)Q 
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In three dimensions, the integral over Q exhibits an ultraviolet divergence arising from the failure of the diffusion 
equation ([5]) at small length scales |r — r'| < I. This unphysical divergence can be regularized by introducing an upper 
cutoff of integration Q m ax ~ 1 /(■■ 

Although, strictly speaking, Eq. §5§ with Db replaced by £>(ft) can only be justified for kl ^S> 1, the great success of 
self-consistent equations (|5|) and ([6]) is due to the fact that they correctly describe many aspects of wave propagation 
in disordered media all the way down to kl ~ 1 (mobility edge) and even at kl < 1 (Anderson localized regime). 
In a disordered metal, for example, where the quantity of interest is the dynamic conductivity cr(ft) oc -D(ft), these 
equations yield the weak localization effect cr(0) oc 1 — const/(fc£) 2 , the low-frequency behavior of conductivity at 
the mobility edge cr(ft) oc (— ift) 1 / 3 and in the localized (i.e. insulating) phase <r(f2) oc — ift£ 2 15:]. However, not all 
the results obtained in the framework of SC theory are correct. As an example, we mention the critical exponent 
v describing the divergence of localization length £ with kl — 1: £ oc \kt — 1| _I/ - SC theory yields v = 1, whereas 
numerical simulations suggest v rs 1.5 (l6| . Another shortcoming of SC theory is its inapplicability to systems with 
broken time-reversal symmetry. 

The derivation of Eq. (|6|) heavily relies on the translational invariance and cannot be straightforwardly generalized 
to media of finite size, even when the size L of the medium is much larger than A and I. To some extent, Eq. 
§5§ with Db replaced by -Dfft) can still be used to study media of finite size by using a lower cutoff ~ \/ L in the 
integral over Q in Eq. ([6]) [5J. Such an approach can be more or less successful in making qualitative predictions in 
the spirit of the scaling theory of localization [T^, but it becomes insufficient when one is interested in fine details 
of multiple wave scattering close to the mobility edge and in the localized regime: coherent backscattering cone 6j, 
dynamics of short pulses @, [1| , or precise scaling of the transmission coefficient with the size of disordered sample 9| . 
A plausible generalization of SC theory to media of finite size can be obtained by noticing that, by virtue of Eq. (0), 
the Q-integral of Eq. © is formally equal to the "return probability" G(r, r, ft). We can therefore rewrite Eq. ^ as 
1/D(ft) = 1/Db + 6ir/ (k 2 l)C(r, r, ft). Van Tiggelen et al. conjectured 6] that in this new form the self-consistent 
equation for D might hold in a medium of finite size as well. In a medium of finite size, the position dependence of 
C(r, r, ft) naturally gives rise to a position dependence of D: 

1 1 67T 

^ + 7^(r,r,Q). (7) 



D(r,ft) D B k 2 l 

If we then enforce diffusive behavior of the intensity Green's function and insist on the energy conservation, the 
equation for C becomes @ 

(-in - V r ■ D(r, ft) V r ) C(r, r', ft) = S(r - r'). (8) 

Although SC equations ([7]) and {5} appear to be a powerful tool to study Anderson localization in realistic situations 
[H, 0, B @|j they still remain a conjecture and lack microscopic justification. Derivation of these equations from the 
first principles is the main purpose of the present paper. 

III. INTERFERENCE EFFECTS IN FINITE MEDIA 

Formally, the intensity Green's function is given by [23| 

4-7T 

C7(r,r',ft) - — (G(r,r', Wl ))(G*(r,r',a; 2 )) 
c 

47T 
+ 

c 

x r(r 1 ,r 2 ,r 3 ,r 4 ,ft)(G(r 2 ,r', Ci ; 1 ))(G*(r 4 ,r', W2 )), (9) 



■ J dridr 2 dr 3 dr 4 (G(r, r 1} (G* (r, r 3 , w 2 )> 



where T(ri, r 2 , r 3 , r 4 , ft) is the complete vertex function given by a sum of all diagrams connecting scattering paths 
corresponding to G and G*. The first term (G)(G*) on the right-hand side (r.h.s.) of Eq. §§§ will be neglected in the 
following. Indeed, it is exponentially small at large distances |r — r'| 3> I that are of main interest for us here. 

In the regime of weak disorder, defined by kl 3> 1, r(r 4 , r 2 , r 3 , r 4 , ft) = <5(r 4 — r 3 )<5(r 2 — r 4 )r£>(ri, r 2 , ft) with a 
sum of ladder diagrams Bam shown in Fig. [TJa). We denote G given by Eq. with To substituted for T by Cd- 
At large distances |r — r'| ^> I and in the limit of small ft, Cd obeys the diffusion equation ([5]). We also introduce a 
sum of maximally-crossed diagrams Tc(yi, r 2 , ft) shown in Fig. [TJb). If we do not consider the first term on the r.h.s. 
of Fig. HJa) , we can formally obtain Tq from by rotating the bottom propagation line of the diagram of Fig. [TJa) 
by 180° in the plane perpendicular to the plane of the figure. The time-reversal invariance, that we assume to hold 
throughout this paper, implies Tc(rx, r 2 , ft) = Yd(t\, r 2 , ft) if |r 4 — r 2 | exceeds the correlation length of disorder (i.e. 
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FIG. 1: (a) Sum of ladder diagrams ro(ri,r2, 0.) and (b) sum of maximally-crossed diagrams Tc(ri, r2, 0,). Solid and dashed 
lines denote (G) and (G*), respectively. The dotted line symbolizes the correlation function of disorder k 4 (fi(r)fi(r')} given 
by Eq. (5). Crosses denote scattering events. Integrations over positions of all internal scattering events are assumed. In all 
diagrams of this paper, (G) and (G*) should be evaluated at frequencies u}\ — ujq + fi/2 and ui2 = ujo — 0/2, respectively. We 
show this explicitly in the panel (a) of this figure only. 




FIG. 2: The diagram X(r, r ,fi) that we use to introduce interference effects in the calculation of intensity Green's function. 
This diagram is made of a four-point Hikami box H(r, n, r', r2) — detailed in Appendix 1X1 — and of the sum of maximally- 
crossed diagrams rc(ri, T2, fi) shown by wavy lines connecting ri and r2. 



if i"i 7^ r 2 for the white-noise disorder that we consider here) because the first term of Fig. [TJa) can be neglected in 
this case. 

To account for interference effects during propagation, we consider a loop-shaped diagram X(r, r', fl) shown in Fig. 
[2j This diagram is made of a square diagram known as a four-point Hikami box H(r, ri, r', r2) [12, uM an( i °f a sum 
of maximally-crossed diagrams rc(ri, r 2 , SI) that we replace by To{^i, i"2, making use of time-reversal invariance: 



X(r, r',Sl) = J dndr 2 iJ(r, n, r', r 2 )T D (r u r 2 , SI). 



(10) 



Because H is a local object having non-zero value only when all the 4 points r, ri, r' and r 2 are within a distance of 
order I from each other, we can expand in series around r, assuming that its spatial variations are small at the 
scale of I: 



r D (n,r 2 ,ft) ~ |l + (n -r)- V ri + (r 2 -r) • V r2 



i[( ri -r).V ri ] 2 + i[(r 2 -r).V r 



}r D (n,r 2 ,r2) | ri= 



(ii) 



We will truncate this expansion to the first order and use the reciprocity principle Y n(?\,ti,Si) = Td{^2, n, Si) that 
allows us to rewrite Eq. (jlip as 



r D (ri,r 2 ,fi) 



l + -(ri+r a -2r)- V, 



T D (r,T,Q). 



(12) 
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Substituting this into Eq. ([TO]) we obtain 



J ff(r,r') + -H / (r,r')-V I 



r D (r,r,fl), 



with 



and 



X(r,r',fi) = 

H(r,r') = J dr 1 dr 2 H(r, n, r', r 2 ) 
H/(r,r') = J dridr 2 (ri+r 2 -2r)iJ(r,n,r',r 2 ) 



(13) 



(14) 



(15) 



The first term on the r.h.s. of Eq. (|13p is the "usual" term arising in the infinite medium as well [l4T |. The second 
term on the r.h.s. is non-zero only in a finite medium because Td(t, r, Q) is independent of r in the infinite medium. 
It will be seen from the following that this term is of fundamental importance for the derivation of self-consistent 
equations of localization in a finite medium. 
A calculation detailed in Appendix [A] gives 



H / (r,r') = -(r-r')^(r,r') 
Substituting Eq. (fl~6)l into Eq. (fl~3)l we obtain 

1 



X(r,r',n) = H(r,r') 



I --(r-r')-V r 



r D (r,r,n). 



(16) 



(17) 



For convenience of calculations, we introduce the difference variable Ar = r — r', such that a given function f o f 
r and r' becomes a function / of r and Ar. In particular, H(r,r') becomes H(Ar) and does not depend on r 
Using the new set of variables r and Ar, we have Tjj(r, r, Q) = T^r, Ar = 0, f2). Equation (jTTJ) becomes 



X{r, Ar,fi) = iT(Ar) 



l-I A r-V r 



r D (r,o,n). 



(18) 



We now take the Fourier transform of Eq. I|18p with respect to Ar and consider the limit q — > 0. Because the Fourier 
transform -ff(q) of if(Ar) is equal to Dgi^q 2 /8irck 2 in this limit [IH, [24j], we obtain 



X(r,q,r») 



-£Di 
8-rrck 2 



[(«q) 2 + (»q) • v r ] f D (r,o,n). 



(19) 



An approximate expression for A(r, Ar, f2) can then be obtained by the inverse Fourier transform of Eq. (fl"9]) with 
respect to q (see Appendix |B|) : 



X(r,Ar,Q) 



-£ 4 De 
8nck 2 



{[A Ar 5(Ar)] + [V Ar S(Ar)} ■ V r }f D (r,0,fi). 



(20) 



Because VAr^(Ar) = V r <5(r — r') and AA r <5(Ar) = A r 5(r — r'), Eq. (|2l7|) can be rewritten in terms of the original 
variables r and r' as 



X(r,r',n) 



~l A D B 
8irck 2 



V r - [r o (r,r,0)V r ]5(r-r') 



(21) 



IV. DERIVATION OF SELF-CONSISTENT EQUATIONS 

We will now use the diagram X of Fig. [5] analyzed in the previous section to include interference effects in the 
calculation of intensity Green's function C(r, r',Q). To this end, we insert the "interference loop" X in the sum of 
ladder diagrams for Cd and account for the possibility of having multiple consecutive interference loops. This leads 
to an infinite series of diagrams shown in Fig. [3] This series can be written analytically as 

Attc f 

C(r,r',Q) = C J) (r,r',!l) + -jC J) (r,r 1 ,!!)I(r ll r 2 ,!l)C J) (r 2 ,r',!!)rfMr 2 
J C D (r, r 1 ,n)X(r 1 ,r 2 ,n)C D (r 2 , r 3 , O) 



/ 47TC 

x X(r 3 ,r i ,n)CD(r4 : ,r\D,)dridr2dr 3 dr i + 



(22) 



FIG. 3: Diagrammatic representation of an infinite series of diagrams contributing to the intensity Green's function. The first 
term is the sum of ladder diagrams. The second term is the sum of ladder diagrams with a single interference loop denoted by 
wavy lines and equal to an infinite sum of maximally-crossed diagrams. Next terms contain 2, 3, etc. consecutive interference 
loops. The ladder and the maximally-crossed diagrams are joined together by a Hikami box detailed in Appendix [X] The 
analytic representation of this diagrammatic series is given by Eq. (|22[1 . 



We now apply the operator — ifl — Z?sA r to Eq. (f2"2")) and use Eq. ([5]) for Cd and Eq. (|2~Tj) for X(r, r'). This yields 
(see the detailed calculation in Appendix [C|l : 

[-tfi - D B A T ] C(r, r', ft) = S(r ~ r') - V r • [T D (r, r, fl) V r C(r, r', ft)] , (23) 



or 



-HI - V,. [D B - ^r D (r,r,0) j V, 



C(r,r',fl)=S(r-r'). (24) 



As we demonstrate in Appendix iDl Yd is proportional to Cd'- rrj(r, r, fl) — (4:1tc/£ 2 )Cd(j, r, ft). This allows us to 
define a renormalized, position-dependent diffusion coefficient 

One 

D(r, ft)=D B - —D B C D {r, r, fl). (25) 

and rewrite Eq. (f2~4"|) as 

{-if} - V r • D(r, fl) V r ] C(r, r', fi) = 5(r - r'). (26) 

The last step consists in applying the self-consistency principle [5[. This can be done by using D(r, ft) instead of 
D B when calculating the second term on the r.h.s. of Eq. (T2"o]) . Diagrammatically, this procedure is equivalent to 
inserting "secondary loops" in the loops shown by wavy lines in Fig. [3] and then inserting the same loops in these 
secondary loops, etc., thus obtaining a sum of diagrams with an infinite sequence of loops inserted one inside the 
other. Physically, this simply means that the same, self-consistent diffusion coefficient D(r, fl) should be used when 
we calculate the intensity Green's function C and the sum of maximally-crossed diagrams Yc- More specifically, we 
have to perform the following replacements: 

1. We replace D B by D in H(r,r') in Eq. (T7|), or equivalently in -ff(q), such that Db is replaced by D in the 
second term on the r.h.s. of Eq. 



2. We replace Db by D in Yd in Eq. (fTT|) , which amounts to replace Cd by C in the second term on the r.h.s. of 

Eq. jjg. 

Equation (25]) then becomes D(r,fl) = D B - {2nc/k 2 )D(r,fl)C(r,r,fl) or 

Dik» = ik + m c{r > r > n) - (27) 

This completes the derivation of self-consistent equations of localization — Eqs. (|2"6")l and (f2"T)) — in a medium of finite 
size. 

The solution of the diffusion equation (T2"6")l in three dimensions diverges when r' — > r: C(r, r', ft) cx l/\v — r'\. This 
unphysical divergence poses potential problems in Eq. (|2T[) that contains C(r, r, ft). One possibility to regularize 
this divergence is to represent C(r, r', ft) as a Fourier transform of C(r, q, ft), where q is a variable conjugated to 
Ar = r — r'. and then cut off the integration over q at some q m ax ~ l/£- The exact proportionality constant between 
q m ax and l/l will determine the exact position of the mobility edge k£ ~ 1. ft is also possible to cut off only the 
integration over qj_ = {q x ,Cly), leaving the integration over q z unrestricted. Such a two-dimensional cutoff is easier to 
implement for the particular geometry of a disordered slab perpendicular to the z axis [H, Q . As could be expected, the 
main qualitative features of final results are largely insensitive to the details of the large-g cutoff, although quantitative 
details can vary slightly. 
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V. ENERGY CONSERVATION AND BOUNDARY CONDITIONS 



It is important to note that although we have obtained Eq. (|26[) by summing only the diagrams of certain type 
and neglecting many other diagrams, this equation satisfies the conservation of energy exactly. Indeed, let us take its 
inverse Fourier transform with respect to ft: 

dC{V d r t '> t] - J g V r • D(r, ft) V r C(r, r', Ct)e~ iat = S(r - r')5(t). (28) 

The flux of energy is given by Fick's law: J(r, r', t) = - J dfl/ (2ir)D(r, ft) V r C(r, r', fl)e~ mt . By integrating Eq. {28"]) 
over a control volume V contained inside the disordered medium and enclosed by a surface S, we obtain 

[ dC ^: r ''^ dr = - ( V r -3(r,r\t)dr + S(t) [ 5(r - r'W. (29) 
Jy at J v J v 

We now apply the Gauss-Ostrogradsky theorem to the first term on the r.h.s. of Eq. (|29|) and assume that the source 
point r' is contained inside V: 



d_ 

dt 



■ [ C(r, r', t)dr = - I J(r, r', i) ■ dS + <5(t). (30) 

Here dS is a vector normal to the surface element dS and directed outwards the volume V. 

Equation (|30[) is a conservation equation. It states that the variation of wave energy in the volume V is given by 
a balance of energy emitted by the source (the second term on the r.h.s.) and energy leaving the volume through its 
surface S (the first term on the r.h.s.). 

Although inside a disordered medium the energy flux J(r, r', t) can have arbitrary magnitude and direction consistent 
with the diffusion equation (|2"rj|) and Fick's law, additional factors come into play at the surface of the medium. More 
specifically, for an open disordered medium of convex shape surrounded by the free space, no energy flux enters 
the medium from outside, provided that all sources are located inside the medium. This simple principle allows a 
derivation of boundary conditions for the intensity Green's function at the surface of disordered medium. Following 
Zhu et al. [221 ] . we consider a disordered medium occupying the half-space z > 0. At a given point r inside the 
medium, the Fourier component of intensity 7(u,r,r',f2) propagating in the direction of a unit vector u, can be 
represented as 0, HI] 

I(u,r,r',Q) = C(r,r',n) + -J(r,r',ft)-u 

c 

= C(v,r',n)--D(r,n)V r C(r,r , ,n)-u, (31) 
c 

where Fick's law was used to obtain the second line. The total flux of wave energy crossing some plane z — const at 
point r in the positive direction of axis z is 

,2tt ,tt/2 

J+(r,r',fi) = — / dd> I d0sin0 u z I(u,r,r', ft), (32) 
4tt Jo Jo 

where u z = cos 9 is the z component of u. We then substitute Eq. (|3ip into Eq. (|32p and perform integrations over 
9 and (f>. This yields 

C(v,v',n)c D(r,Sl) dC(r,r',n) 
J + (r,r,0) = . (33) 

By requiring J+(r, r', it) — - at the surface z = of the medium, we obtain the following boundary condition: 

2 rv m, dC(r,r',rt) 



C(r,r',n)\ z=0 --D(r,n)\ z=0 



c dz 



= 0. (34) 



2 = 



For a medium of more complex but still convex shape, the above derivation can be repeated locally in the vicinity 
of each point of the medium surface S, assumed to be locally flat. This yields 

2 D(r CI) 

C(r, r', CI) - -t^-^- (n(r) ■ V) C(r, r', CI) - 0, (35) 
o L>b 

where n(r) is a unit inward normal to the surface S at the point r £ S. This equation is the boundary condition 
for the intensity Green's function at an open boundary. It can be generalized to include internal reflections of waves 
at the boundary by replacing 21/3 by a larger "extrapolation length" zq in front of the second term on its l.h.s., in 
complete analogy with Ref . |22j. 
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FIG. 4: Hikami box H(r, n, r', r2). Diagrammatic notation is the same as in Fig. \T\ 



VI. CONCLUSION 

In this paper we derived the self-consistent (SC) equations of Anderson localization — Eqs. ([26]) and (|27|) — starting 
from the first principles. Mathematically, this was achieved by dressing the ladder propagator with "interference loops" 
made of maximally-crossed diagrams. Each loop was inserted into the ladder with the help of a Hikami-box diagram. 
The SC equations were then obtained by applying the self-consistency principle. 

The essential difference of our derivation compared to the derivation of SC equations in the infinite medium is 
the position dependence of the sum of ladder diagrams r£>(r,r',f2) with coinciding end points r = r'. This position 
dependence leads to the appearance of an additional term, proportional to V r T D (r, r, f2), in a series expansion of 
Td(ti,T2, f2) around an arbitrary point r. As a consequence, we have to keep an additional term in the expression of 
Hikami box employed to connect ladder and maximally-crossed diagrams in our approach. It is this term that finally 
allows us to derive SC equations of localization in a medium of finite size. 

Although the condition M> 1 was explicitly used to derive Eqs. I|26p and (|27[) . one can still hope that, similarly 
to SC equations in the infinite medium, they could yield reasonable results in the vicinity of mobility edge and in the 
localized regime. According to Refs. @, 0, HI, S this seems indeed be the case. However, one should understand that 
even though the general form of these equations might be largely universal in both diffuse and localized regimes, the 
numerical prefactor 6ir/k 2 £ in front of the second term in the SC equation for D(r, ft), Eq. (|2~T|) , should not be taken 
too seriously because it originates from the calculation of complicated diagrams that was carried out in the limit 
kl 3> 1 only (see Appendix [A"]l . When the result is extrapolated to kl < 1, this prefactor could vary and, in general, 
its dependence on kl is likely to be more complex than just l/(ki) 2 . In addition, the SC theory neglects interference 
processes insensitive to the breakdown of time- reversal invariance by, e.g., a strong magnetic field. The inclusion of 
such processes in the theoretical description would at least change the prefactor in Eq. (f27|) . In Refs. 0, H, H, for 
example, a larger prefactor was used in Eq. (|27p to study the vicinity of the localization transition. This was justified 
by a comparison of some of the final results with those of the supersymmetric a- model (25[. Such a comparison 
indicates that the prefactor 6ir/k 2 £ in Eq. (f2"T| have to be multiplied by 2 to obtain an exact correspondence between 
the two theoretical approaches 0] ■ 

Finally, SC theory of localization is a very convenient tool for description of realistic experimental situations, like 
the recent experiments on Anderson localization of light [26], microwaves (27|, ultrasound (28[, and matter waves 
[IH [l9j . It can be adapted to almost any detail of a particular experiment (short pulses or focused beams, internal 
reflections on the sample surface, complex shapes or inhomogeneous scatterer density profiles of disordered samples, 
etc.). This gives SC theory a serious advantage as compared to other theories of Anderson localization. 

Note. After this paper was submitted for publication, we became aware of the work of C. Tian (29| who justifies 
the concept of the position-dependent diffusion coefficient using methods of supersymmetric field theory. 
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APPENDIX A 

In this appendix we present a demonstration of Eq. (jTHJ) . Consider Eq. (fT5|) , where the Hikami box H(r, ri, r', T2) is 
shown Fig.[D This diagram is a sum of three contributions: H^ A '(r, Ti, r', T2), H^ b ' (r, ri, r', v?) and H^°'(r, ri, r', Tz)- 



9 



We hence have to perform three integrals. The second one, for example, is 

Hj, B) (r,r') = y'dr 1 dr 2 (r 1 +r 2 -2r)i?( s )(r,r 1 ,r / ,r 2 ). (Al) 

The two other integrals H^ 4 -* and are defined similarly. In the following we focus on the calculation of H^ 8 -*, 

the calculation being similar for Hy 4 " 1 and H^ 7 -* . Eq. (|A1|) can be rewritten as 

Hf } (r,r') = ^ |dr 1 dr 2 dr 3 (r 1 +r 2 -2r)(G(r,r 3 ))<G(r3,r 1 ))<G*(r 1 ,r')) 

x (G*(r,r 2 ))(G(r 2 ,r 3 ))(G(r 3 ,r')). (A2) 

In this Appendix, we neglect the difference in frequencies lo± and w 2 in the arguments of (G) and (G*), respectively, 
and set uj\ = oj 2 = loq for all amplitude Green's functions. This is justified as far as slow dynamics (57 = u>i — oj 2 <C u>o, 
c/tj is concerned. To lighten the notation, we omit the frequency argument of (G). 
By replacing the Green's functions in Eq. (|A2[) by their Fourier transforms, we obtain 



Hf^r,!-') - J dr 1 dr 2 dr 3 dk 1 ...dk 6 (G(k 1 ))(G(k 2 ))(G*(k 3 )) 

x (G*(k 4 ))<G(k 5 ))(G(k 6 ))(r 1+ r 2 -2r) 

x e -ir3(ki-k2+k5-k 6 ) e -iri(k2-k 3 ) e -ir2(k4-k5) e -ir(-ki-k4) e -ir'(k3+k 6 ) 

= K 1 (r,r')+K 2 (r,r , ) + K(r,r / ), (A3) 

where Ki(r, r') is the part of Eq. (|A3|) with the integrand proportional to rj., K 2 (r, r') is the part with the integrand 
proportional to r 2 , and K(r, r') is the one with the integrand proportional to — 2r. Let us first consider Ki(r,r'). 
In this term, the integrals over r 2 and r 3 give respectively (27r) 3 (5(k 4 — k§) and (27r) 3 <5(ki — k 2 + k5 — ke), and the 
integral over ri gives — i(27r) 3 V ka (5(k 3 — k 2 ). We have then 

— 47T? f 

Ki(r,r') = j^J dk 1 dk 2 rfk 4 (G(k 1 ))(G(k 2 ))(G*(k4))<G(k4))(G(k 1 -k 2 +k 4 )) 

x e- ir( - kl - kl) ydk 3 (V k3 5(k 3 -k 2 ))(G*(k 3 ))e- tr ' (k3+kl - k2+kl) . (A4) 
The integral over k 3 is equal to -((Vk 2 (G*(k 2 ))) - ir'(G*(k 2 ))) e -""'( kl + k4 ) and hence 

4-7T7 r 

Ki(r,r') = J dk 1 dk 2 dk 4 (G(k 1 ))(G(k 2 ))(G*(k 4 )) 

x (G(k 4 ))(V k2 (G*(k 2 )))(G(k 1 -k 2 + k 4 ))e^- r ')( kl+k *) 
+ ^9 / dk 1 dk 2 dk 4 (G(k 1 ))(G(k 2 ))(G*(k 4 )) 

x (G(k 4 ))(G*(k 2 ))(G(k 1 -k 2 +k 4 ))e^- r '» kl+k ^. (A5) 
The second term on the r.h.s. is nothing else than r'H^^r, r'), where 

H^(r,r')= J dridr 2 ff( s )(r,n,r',r 2 ). (A6) 

In the first term on the r.h.s. of Eq. (|A5j) we change the variables k 4 — > k, k 2 — > k', and k 4 + k 4 — ► q. Equation 
151) becomes 



47T7 r 

Kl(r ' r ' } = J{^f I dk rf k '^(G(k))<G(k'))(G*(q-k)) 

x (G(q-k))(V k -(G*(k , )))(G(q-k'))e tq(r - r ' ) +r'ij( B )(r,r'). (A7) 
In the limit of small q, we have 1/(2tt) 3 / dk(G(k)) (G*(q - k)) (G(q - k)) = -i£ 2 (l - q 2 ( 2 /3)/%nk [3 and 

Kl (r,r') = J^r [ dqe^-^- ( 1 - 

V ' ; (2tt) 3 J H 2k V 3 / 

x j±- J dk'(G(k'))(G(q-k'))(V k ,(G*(k')»+r'ff( B )(r,r'). (A8) 
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A similar calculation gives 

K ^'> = -<^/^ <,< "' , s( 1 - i ?) 

* (^/' ik '< G ( k '»< G ( c l- k '»(V k .(G-(k')))+r//l'"(r,r'). (A9) 

Besides, it follows straightforwardly from Eq. (|A2j) that K(r, r') = — 2rij( B )(r, r'). Combined with Eqs. (|A8jl and 
"T9]), this yields 



Hf } (r,r') = Ki(r,r')+K 2 (r,r')+K(r,r') 

= -(r-r')ii" (B) (r,r'). (A10) 

The calculation of H (A) and H (c) follows the same lines. We obtain 

H ( f A) (r,r') = J 'rfr 1 rfr 2 (r 1 +r 2 -2r)^ A )(r,r 1 ,r',r 2 ) 

= -(r-r')# (A) (r,r') (All) 

and 

Hf } (r,r') = y'dr 1 dr 2 (r 1 +r 2 -2r)ff^(r,r 1 ,r',r 2 ) 

= -(r-r')-ff (C) (r,r')- (A12) 
Combining Eqs. fArOf . (fATT|) . and (|AT2l) we find 

H/Cr,^) = HW(r,r') + Hf'(r,r')+Hf'(r,r') 

= -(r-r')^(r,r'), (A13) 

which is Eq. (I16|) of the main text. 

APPENDIX B 

We show here that the inverse Fourier transform of Eq. (|19| with respect to q is given by Eq. (|20p . As a function 
of q, Eq. (fT5|) is a sum of two terms proportional to (iq) 2 and «q, respectively. The inverse Fourier transform of iq is 



dq zq e iq Ar = V Ar 



(27T) 



<iq 



= V Ar <5(Ar). (Bl) 



Similarly, the inverse Fourier transform of (iq) 2 is AA r <5(Ar). This leads directly to Eq. (1201) . 



APPENDIX C 



Here we obtain Eq. (|23[) from the series of Eq. (f2"2"|) . The idea is to apply the operator — ifl — -D_eA r to both sides 
of Eq. ([22]) . The first term on the r.h.s. is transformed into S(r — r') since Co(r, ri, fi) obeys Eq. ([5]), and for each 
of the next terms the first multiplier Cd(t, ^) in the integrands is transformed into <5(r — ri) for the same reason. 
Equation (|2"2")) becomes 

[-iO-D B A r ]C(r,r',{l) = 6{r - r') + — / <5(r - ri)X(n, r 2 , Q)C D (r 2 , r', fi)dndr 2 

+ fe) |^r-r 1 )A(r 1 ,r 2 ,f])^(r 2 ,r 3 ,r!) 

x A(r3,r4,f])C D (r4,r',n)dridr 2 dr 3 dr4 + . . . (CI) 
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Performing integrations over ri we obtain 

4-7TC f 

[-iSl-D B A r ]C(r,r',Sl) = 6(r - r') + — / X(r, r 2 , Sl)C D {v 2l r', Sl)dr 2 

+ (7?) / X(v,v 2 ,Sl)C D (v 2 ,r 3 ,Sl) 

x X(r 3 ,r4,n)C D (r4,r^)dr 2 dr 3 dr 4 + ... (C2) 
Now let us perform integrations over r 2 . We have to calculate an integral 

4irc 



r 

X(r,r 2 ,Sl)C D (r 2 ,r 3 ,Sl)dr 2 , (C3) 
where r3 = r' for the first integral on the r.h.s. of Eq. (|CJ2|) . To this end, we use Eq. f2"T)) for X(r, r 2 , SI) and obtain 

" ' ~ DB J dr 2 [V r • T D {r, r, Sl)V r 6(r - r 2 )] C D (r 2 , r 3 , fi) 

i;2 -V r -[r D (r,r,fi)V r ]C B (r,r 3 ,r!), (C4) 



2fc 2 
2fc 2 



where we integrated by parts. 

Replacing each integration over r 2 by this result in Eq. (|C2[) we obtain 



[-ifi - D B A r ] C(r, r', 0) = 5(r-r') 

: -V r -r D (r,r,n)V r \c D (r,r',n) 



2k 2 



+ — / ^(r,r 3 ,ri)X(r3,r4,r!)C I) (r4,r',fi)dr3dr4 

x X(r 5 ,r 6 ,Sl)C D (r 6 ,r',Sl)dr 3 dr 4 dr 5 dr e + ... (C5) 

The infinite series in square brackets is nothing else than the intensity Green's function C(r, r', SI) as given by Eq. 
((22)1 . Thus, Eq. (C5j| leads straightforwardly to Eq. ((25)1 . 

APPENDIX D 

We prove here that the proportionality between Cu an( i known in the infinite medium [l4| holds in a finite 
medium as well. Calculations being similar to those of Appendix (X) we only give the main ingredients of the proof. 
According to Eq. ©, at |r — r'| 3> t, C l £>(r, r', SI) is given by 

47T f 

C D (r,r',Sl) = — I dT.dr^Gir^^iG^r^^rnin^^tyiG^^'ViG^r^r')). (Dl) 

Similarly to Appendix (X) we omit frequency arguments of amplitude Green's functions and set all of them equal to 
ojq. Because (G(r, ri)) is exponentially small for |r — ri| > £, the main contribution to the integral comes from |r — 1*1 1 , 
|r' — r 2 | < I. This authorizes us to expand Td(yi, r 2 , SI) in series around (r, r'). This expansion has to be truncated 
to the same first order in | r — r 1 1 and |r' — r 2 | as the expansion of Eq. (fT2"]l : 

r D (n, r 2 , Si) ~ T D (r, r', Si) + (n + r 2 - r - r') • V r T D (r, r', SI). (D2) 

We then substitute Eq. (|D2p into Eq. (|D1[) . The integral proportional to Ti)(r,r' ,Sl) is the usual result obtained in 
the infinite medium. This integral equals ^ 2 /(47rc)r£>(ri, r 2 , SI) [14(. We hence obtain 

C D (r,r',Sl) = ^—T D (r,r',Sl) (D3) 

47TC 



47T 
+ 

c 



J dv 1 dv 2 {v 1 +v 2 - v - r') |(G(r,n))| 2 |(G(r 2 ,r'))f 



V t T d (t,t',SI). 
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The integral on the r.h.s. of Eq. (|D3|1 can be calculated exactly in the same way as Hf , or in Appendix [Al 
and it is easy to see that this integral is zero. 
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